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ABSTRACT 

We determine the characteristics of the 7 mm to 20 cm wavelength radio variability in Sgr A* on 
time scales from days to three decades. The amplitude of the intensity modulation is between 30 and 
39% at all wavelengths. Analysis of uniformly sampled data with proper accounting of the sampling 
errors associated with the lightcurves shows that Sgr A* exhibits no 57- or 106-day quasi-periodic 
oscillations, contrary to previous claims. The cause of the variability is investigated by examining 
a number of plausible scintillation models, enabling those variations which could be attributed to 
interstellar scintillation to be isolated from those that must be intrinsic to the source. Thin-screen 
scattering models do not account for the variability amplitude on most time scales. However, models 
in which the scattering region is extended out to a radius of 50-500 pc from the Galactic Genter account 
well for the broad characteristics of the variability on > 4-day time scales. The ~ 10% variability on 
< 4-day time scales at 0.7 — 3 cm appears to be intrinsic to the source. The degree of scintillation 
variability expected at millimeter wavelengths depends sensitively on the intrinsic source size; the 
variations, if due to scintillation, would require an intrinsic source size smaller than that expected. 
Subject headings: galaxies: active — Galaxy: center — scattering 


1. INTRODUCTION 

The compact radio source associated with the black hole at the Galactic Center, Sgr A*, is known to vary at 
millimeter and centimeter wavelengths on time scales from hours to years (Brown & Lo 1982; Zhao et al. 1992; Bower 
et al. 2002; Herrnstein et al. 2004). The origin of these variations remains unclear, with strong arguments for both 
extrinsic and intrinsic mechanisms having been advanced (e.g. Zhao et al. 1989 compared with Zhao et al. 2001). 

Interstellar scintillation is the primary mechanism which may cause any extrinsic variability. The same plasma that 
is responsible for the scatter-broadening of Sgr A* at millimeter and centimeter wavelengths (e.g. Lo et al. 1998, 
Bower et al. 2004) is also expected to cause the source to exhibit refractive intensity variations. It has been argued 
that much of the monthly to yearly variability in Sgr A* at wavelengths longer than 6 cm can be explained in terms of 
refractive interstellar scintillation provided that scattering material moves across our line of sight at a relatively high 
speed of ^ 1000— 2000kms“^ (Zhao et al. 1989). However, the variability amplitude is not so easily accounted for: 
a scattering medium modeled as a single thin screen underpredicts the observed variability amplitude (ibid.) while 
extended medium models, which are in principle capable of explaining higher refractive modulation amplitudes for the 
same degree of scatter broadening, have not been investigated in the context of the Galactic Center. 

Recent interpretations favor an intrinsic origin for much of the centimeter wavelength variability. These center 
around claims of 106-day quasi-periodic variations at wavelengths shorter than ~ 3 cm (Zhao et al. 2001) and of 57- 
day quasi-periodic behavior at 2.3 GHz (Falcke 1999). The oscillations possess only a modest spectral purity, with the 
highest purity vj /S.v = 2.2±0.3 reported at 1.3cm. Zhao et al. (2001) discuss the origin of these oscillations in terms of 
periodic flares from a jet nozzle or an instability in the accretion disk triggering, for example, quasi-periodic production 
of convection bubbles. It is widely supposed that the oscillations must reflect a process intrinsic to Sgr A* itself because 
scintillation is incapable of producing such regular oscillations. Yet observations of certain intra-day variable quasars, 
whose variations are proven to be scintillation-induced, invalidate this argument because their fluctuations often exhibit 
even higher degrees of spectral purity (Kedziora-Chudczer et al. 1997; Rickett, Kedziora-Ghudczer & Jauncey 2002). 

Nonetheless there is little dispute that at least some of the variability is intrinsic. Detections of flares at millimeter, 
IR and X-ray wavelengths (Wright & Backer 1993; Tsuboi, Miyazaki & Tsutsumi 1999, Eckart et al. 2004; Baganoff 
et al. 2001) conclusively demonstrate that the source is intrinsically variable. A possible connection between X-ray 
flaring and unusually large flux density excursions at 7 mm is also reported (Zhao et al. 2004). However, it is difficult 
to ascertain how much variability observed at centimeter wavelengths could be attributed to flaring since neither the 
duty cycle nor the energy distribution of mm or X-ray flares is well-constrained, much less the physical connection 
between centimeter and mm or X-ray behavior. 

Despite the many recent observational results concerning the properties of Sgr A*’s variability, a dearth of corre¬ 
sponding theoretical efforts has failed to place these results in context, leaving us none the wiser as to their cause. 
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For instance, while it is acknowledged that scintillation variability is likely to be important at centimeter wavelengths, 
no variations have been specifically attributed to it, and no realistic modeling has been applied to investigate what 
contribution it could conceivably make. This paper aims to redress the balance by investigating two outstanding issues: 
(i) what exactly does a model of Sgr A*’s variability need to explain and (ii) can one deduce which variations must 
be intrinsic to the source by eliminating the variations that can be explained by interstellar scintillation? The next 
section of this paper is devoted to the former question, including a critical examination of the ~ 100 day quasi-periodic 
oscillations reported in Sgr A* (Zhao et al. 2001), while lEI addresses the latter question. We compare the models to 
the observations in 0 and summarize our findings and briefly detail their implications in 0 

2. DATA ANALYSIS 

Sgr A* has been the subject of numerous VLA monitoring campaigns since 1975. The resulting data are published 
in Zhao et al. (1992, 2001) and Herrnstein et al. (2004). The latter lists the results of a three-year effort to measure 
weekly variations at 7 mm, 1.3 cm and 2 cm. We combine all these data to form lightcurves at 7 mm, 1.3 cm, 2 cm, 
3 cm, 6 cm and 20 cm in an attempt to quantify the variability of Sgr A* on time scales of a few days to decades. 

Additional daily Green Bank Interferometer (GBI) monitoring at 2.3 and 8.3 GHz (Falcke 1999) quantifies variations 
in Sgr A* on shorter time scales. We reanalyze these data here, but do not incorporate them with the VLA flux 
density measurements because the GBI is highly susceptible to confusion in the Galactic Center region. The GBI is a 
two-element interferometer whose 2400 m spacing is insufficient to resolve out much of the extended emission near Sgr 
A* which, if not properly accounted for, can cause hour-angle dependent variations in the measured flux density of Sgr 
A*. Falcke (1999) attempted to correct for hour-angle dependent gain variations, and to eliminate the contribution 
of confusion by comparing GBI flux density measurements with available contemporaneous VLA measurements. But 
only by comparing the visibilities to a complete synthesis image of the crowded Galactic Center region can one be 
confident in removing the effect of confusion. 

The lightcurves from the combined datasets are shown in Figure ^ Various parts of the lightcurves have been 
published elsewhere, and their main purpose here is to illustrate exactly which data are used here. 

Throughout this paper we adopt the structure function as the chief measure of source variability because it is ideally 
suited to the interpretation of data which are highly irregularly sampled in time. The intensity structure function, 
Dj{t) = {[I{t' + t) — /(t')]^), is a simple statistic which characterizes the variance between measurements separated 
by a time interval r. Since we wish to be confident that we interpret only those features of the variability that are 
statistically significant, this requires a rigorous assessment of the errors associated with our measure of variability, and 
a statistic simple enough that the errors are readily calculable (see H2.1l belowL The power spectrum, which is related 
to the structure function by a Fourier transform, is a more elegant measure of variability, but we do not employ it 
here because the irregular time sampling of our datasets complicates the error analysis and thus the interpretation of 
the statistic. 

In computing a single structure function to characterize variability over the entire lightcurve we make the implicit 
assumption that the variability statistics are wide-sense stationary, which is to say that the statistical properties 
of the variations themselves do not vary with time. This approach is not strictly valid, for instance, if the source 
undergoes various “phases” of variability in which the presence of fast time scale variations is modulated by some 
underlying long-term process; X-ray binaries, whose behavior is characterized by infrequent outbursts, represent an 
obvious counterexample. Although the possibility that variability in Sgr A* changes character with time cannot be 
discounted, there is no strong evidence to support the notion. 

2.1. Assignment of Errors 

The correct determination of errors associated with any measure of the variability is crucial in assessing the signif¬ 
icance of variability time scales or even periodicities in the data, particularly when they are irregularly sampled and 
when the time scales under consideration are comparable to the entire data length. It is also crucial when comparing 
the observed variations to theoretical models, as most theories only predict ensemble-average quantities (i.e. they 
predict the climate, not the weather). 

The largest contribution to the error arises because our observations only sample the stochastic fluctuations over 
a finite duration. A simple argument would suggest, for instance, that in a dataset spanning 1000 days a process 
that operates on a time scale of 200 days contains at most five independent measurements of the variations caused by 
this process. Formal arguments show that the error in the structure function at delay r from an observation of total 
duration T is (Jenkins and Watts 1968, see also Appendix B of Rickett, Coles & Markannen 2000) 

4 pT-t 

var[D{T)] = _ |^|^2 [j{\r - t\) + -f{\r + t\) - 2j{r)f (T - T - r)dr, (1) 

where the function 7 (t) is the ensemble-average autocovariance of the intensity fluctuations at wavelength A. This 
function is unknown, so we assume that the measured autocovariance is a reasonable representation of its corresponding 
ensemble-average counterpart. 

2.2. Characteristics of the variability 

The structure functions derived from the lightcurves are displayed in Figures Inland 13 Figure [3 is the main result 
of this section, but the plots in Figure|3 shown on linear scales, allow closer scrutiny of the variability characteristics 
on time scales shorter than 200 days. 
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Fig. 1.— Lightcurves of the variations in Sgr A* from 7 mm to 20 cm. Flux density errors were not included in the GBI data reduced by 
Falcke (1999), but the spread of the data points at each epoch is a reasonable indicator of the uncertainty. 


2.2.1. Intra- to inter-day fluctuations 

Sgr A* is reported to exhibit variations on scales down to less than one day (Brown & Lo 1982). It is possible to 
determine whether the present datasets show evidence for intra-day (< 4-day) flickering by examining the behavior of 
the structure functions at small time lags. Short time scale flickering is present if the value of the structure function in 
the smallest time bin differs significantly from zero. A proper assessment of the presence of flickering depends crucially 
on the correct determination of the errors associated with the flux density measurements. We have subtracted the 
contribution of measurement errors from each structure function using the errors quoted in the papers from which the 
observations are derived^. However, if the errors are under-estimated the structure function is biased towards high 
values, leading to an over-estimate of the amplitude of short time scale flickering. Conversely, an over-estimate of 
the errors biases the structure function to lower values, the most clear manifestation of which is to cause dips in the 
structure function below zero, which is clearly unphysical. 

The structure functions in Fig.|31 suggest that Sgr A* undergoes appreciable flickering on < 10-day time scales at 
wavelengths from 7 mm to 3 cm. In all cases the root-mean-square fluctuations are approximately 10% of the mean flux 
density. Inspection of more finely-binned structure functions reveals that Sgr A* exhibits 10%, 6% and 8% variability 
at 7, 13 and 20 mm respectively on < 4-day time scales. However, we caution that our estimate of the flickering 
amplitude depends critically on a correct assessment of the errors associated with the observations. Such an estimate 
is called into question for the 20 cm lightcurves, for instance, where the estimated contribution of measurement errors 
leads to negative structure function values at certain lags. 

The amplitude of the intra-day variability must be subtracted when comparing the observed structure functions to 
models which only apply to variations on longer time scales. However, the amplitude of any intra-day variation is 

^ The means of removing measurement errors is obvious when one writes the structure function in terms of the autocovariance: D(t) = 
2[C(0) — C(t)], where C(t) = r)A/(t)}. Measurement noise is assumed to be uncorrelated between samples and independent 

of the true intensity fluctuations, so it only makes an (additive) contribution to C(0). It is removed by subtracting the variance of the 
measurement errors from C(0). 
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Fig. 2.— Structure functions of tlie variations in Sgr A* from wavelengths of 7mm to 20 cm. The continuous and dashed error lines 
indicate the 1- and 2 -(t confidence limits of the measured structure functions due to the finite duration of the lightcurves (see il2dj. The 
ticks associated with individual points display the error caused by the finite number of observations that contributed to the measurement 
at that particular time lag. 


small compared to total amplitude of the intensity variations, and this correction is small compared to the uncertainty 
in the total variability amplitude on long time scales at most wavelengths. 

2.2.2. Quasi-periodic variations 

Figure El can also be used to assess whether any variations exhibited by Sgr A* are quasi-periodic. Such variability 
is characterized by the presence of oscillations in the structure function. The origin of this behavior is understood by 
noting that the structure function is related to the power spectrum of the lightcurve, P{uj), by a Fourier transform: 
Dj(t) oc J[1 — exp(ia;r)]P(a;)dai. A purely sinusoidal signal in the lightcurve would manifest itself as a sharp peak 
in the power spectrum, and the structure function would exhibit a peak at a time scale corresponding to the period, 
followed by sinusoidal oscillations peaking at multiples of the fundamental period. The amplitude of these oscillations 
would decrease with time lag if the variations are spectrally impure. 
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Fig. 3.— Structure functions of the variations in Sgr A* from wavelengths of 7mm to 20cm. For small lags At < 200 days the errors 
are dominated by the finite number of observations that contributed to the measurement at that particular time lag. The contribution of 
measurement errors is estimated using the errors quoted from the observations, and has been subtracted from each structure function. At 
6cm and 20cm these clearly overestimate the true error, as these structure functions are negative at certain time lags. 
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Fig. 4.— A simple illustration of the effect of quasi-periodic oscillations (top) and flaring (bottom) on the behavior of the structure 
function. The left panels show the sample lightcurve, the middle panel its power spectrum and the right panels the resulting structure 
function. The quasi-periodic oscillations possess a relative bandwidth Azz/z/ of 0.5. 


It is important to distinguish between oscillations in the structure function from spikes which are devoid of accompa¬ 
nying oscillations at longer time scales. A structure function containing sharp, isolated spikes indicates that the power 
spectrum contains quasi-periodic features. This in turn indicates the presence of sharp spikes in the corresponding 
lightcurve. The structure function contains spikes (or sharp dips) when at least two flares are present in the lightcurve. 
For instance, flares at times ti and t 2 each of duration At would give rise to a feature in the structure function at 
a time lag — ti| of width 2At. In the absence of any quasi-periodic or flaring behavior the structure function is 
expected to increase monotonically with time until it saturates at the longest time scale of the variations present in 
the data. Figure 0] illustrates how quasi-periodic oscillations and flares are manifested in structure functions. 

We investigate the form of the structure function on < 150 day time scales, motivated by reports of quasi-periodic 
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variations on time scales of ~ 106 days between 7 mm and 3 cm (Zhao et al. 2001). The structure functions derived 
from the long duration VLA datasets (Fig. OJ are most useful in assessing the data for the presence of any unusual 
features. A simple test for the significance of any features is obtained by fitting a single line through each structure 
function at time lags from At = 5 — 150 days and computing the reduced statistic, as listed in Table 1. The 
departure of the structure function from a line, indicated by a high signifies the presence of peaked features above 
the generally-increasing trend with time. The statistics in Table 1 show that the 7 mm structure function is well- 
fit by a single line, and there is no significant detection of any quasi-periodic variability. However, the reduced 
statistic suggests the presence of significant deviations at all lower wavelengths, as is obvious by inspection of Fig. |3 
The 1.3, 2 and 3 cm structure functions appear to exhibit peaks at lags of At = 30 — 45, 45 — 65 and 45 — 55 days 
respectively. The 2 cm structure function also exhibits a peak at At = 80 — 95 day which appears to coincide with 
a marginally significant peak at a corresponding time lag in the 7 mm structure function, however the coincidence 
does not necessarily increase the significance of the peaks because the errors bars at the two wavelengths are not 
independent if the two lightcurves are partially correlated. In addition, the 6 and 20 cm structure functions possess 
highly significant peaks at At 140 — 170 days; the 1.3 cm structure function appears to exhibit a feature at similar 
time lags, but inspection of the error bars in Fig. [^suggests the detection is of marginal significance. 

We also consider the lightcurves measured with the GBI. Our reanalysis of the GBI dataset reproduces the structure 
functions reported by Falcke (1999) on whose basis a 57-day quasi-periodic variability cycle is claimed at 2.3 GHz. 
These variations were suggested to be quasi-periodic because the structure function subsequently oscillates weakly after 
peaking at ^ 57 days. However, our estimate of the errors associated with these structure functions casts doubt on the 
significance of any claim of quasi-periodic behavior. At both 3.6 and 13 cm even the l-u error troughs are consistent 
with structure functions that increase monotonically and subsequently saturate without undergoing any oscillatory 
behavior. The insignificance of the quasi-periodic behavior is not affected by changes in the temporal binning of the 
structure function. The error troughs reflect only the error incurred by trying to infer the ensemble average behavior 
of the variations from a dataset of finite duration. 

In summary, all of the particular features identified in the VLA structure functions are single, isolated peaks. None 
of these can be attributed to oscillatory behavior. As a further, more sensitive test for quasi-periodic variability we 
present in Fig.^lLomb periodograms for 2 cm, 1.3 cm and 0.7 cm lightcurves based on well-sampled data from 2000 
to 2003. We compare the power spectral density for the data against the 99th percentile expectation of uniform noise 
(dot-dashed curve) and noise with a red spectrum (dashed line). The red spectrum is calculated using 300 Monte 
Garlo simulations that use the sampling function of the data sets. We see clearly that there are no significant periods 
in the data. In particular, there are no peaks in the vicinity of the 106-day period reported by Zhao et al. (2001). The 
spikes in the 2 cm and 1.3 cm data at ^ 0.003 days”^ have a significance of only ^ 90% for the uniform noise case. 
Moreover, these periods have only been sampled ~ 3 times in this data set and they are not apparent in PSDs from 
longer (but less well-sampled) data sets. Against the red noise case, these peaks have minimal significance. Results 
are similar to the uniform noise case when one estimates the PSD significance through Monte Carlo simulations in 
which the lightcurves are generated through reordering of the data. We conclude, then, that there is no evidence for 
periodic or quasi-periodic oscillations in the radio light curves of Sgr A*. 

2.2.3. Essential characteristics of the variability 

The structure functions presented in Fig. |21 are the key observable that any variability theory must reproduce. 
Broadly speaking, each structure function may be characterized by a monotonically increasing portion until it saturates 
at a certain amplitude and time scale. A viable explanation of the variability should explain the shape of the structure 
function over an appreciable range of time lags. Even if the model does not explain every significant peak and wobble 
in the structure function, it is still viable if it explains the (i) amplitude and (ii) time scale at which the structure 
function saturates and (iii) the slope of the structure function over an appreciable range in time lags. In the following 
sections we gauge the success of our model by its ability to reproduce these three characteristics. 

As remarked above, many of the structure functions exhibit spikes. These are caused by large, rapid and isolated flux 
density excursions in the lightcurves. We shall not attempt to explain these features in this paper. We merely remark 
that they likely represent either flaring activity intrinsic to Sgr A* itself, or a manifestation of Extreme Scattering 
Events (Fiedler et al. 1987) which appear to occur when the Earth traverses a caustic surface of certain lens-like 
discrete objects in the interstellar medium (Romani et al. 1987; Walker & Wardle 1998). The prevalence of these 
flares in most of the structure functions indicates that the mechanism responsible operates over a factor of ten range 
in wavelength. 

3. SCINTILLATION VARIABILITY 

In this section we attempt to distinguish between variations intrinsic to Sgr A* itself and those due to refractive 
interstellar scintillation. Several detailed scintillation models that span the range of possible scattering conditions are 
constructed in order to isolate those variations which cannot be explained under any plausible scattering conditions, and 
must therefore be attributed to intrinsic source variability. The distinction between intrinsic and extrinsic variability 
is made on the basis of a scintillation model because the physics of any centimeter-wavelength intrinsic variability in 
Sgr A* is ill-constrained; indeed even the fraction of the radio emission originating in the jet and accretion disk is 
disputed (e.g. Falcke & Markoff 2000; Quataert & Narayan 1999; Yuan, Markoff, Falcke 2002). On the other hand, 
the basic physics of interstellar scintillation is well-understood and makes robust predictions that can be compared 
directly to the observed intensity variations. 
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Fig. 5.— Lomb periodograms of the variations between 2000 and 2003 at 20, 13 and 7mm. The power spectral density for the data 
(solid line) is plotted with the 99th percentile expectation of uniform noise (dot-dashed line) and noise with a red spectrum (dashed line). 
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0.05 
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2 

0.05 
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0.1 
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1.0 ±0.2 
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20 
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1.6 ±0.6 
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TABLE 1 

Characteristics of the structure functions in Figures El The index of the power law in the increasing part of the 

STRUCTURE FUNCTION IS OBTAINED FROM A FIT OF THE FUNCTION A + Bt°‘ TO THE DATA BETWEEN DAYS 5 AND 40 IN THE 
LOGARITHMICALLY-BINNED STRUCTURE FUNCTIONS. ThE CONSTANTS B AND a ARE ALLOWED TO VARY, BUT THE OFFSET A IS CHOSEN TO 
REMOVE ANY CONTRIBUTION THAT REFLECTS INTRA-DAY VARIATIONS WHOSE CONTRIBUTION TO THE STRUCTURE FUNCTION DOES NOT 
APPEAR TO BE ASSOCIATED WITH THE VARIATIONS ON TIME SCALES LONGER THAN 5 DAYS (i.E. INTRA-DAY VARIABILITY); THIS PARAMETER IS 
NONZERO AT 7, 13 AND 20 MM ONLY (SEE THE DISCUSSION IN il2.2.1l . ThE LAST COLUMN LISTS REDUCED VALUES RESULTING FROM 
ATTEMPTS TO FIT EACH STRUCTURE FUNCTION AT LAGS BETWEEN 5 AND 150 DAYS BY A STRAIGHT LINE. 


The shortcoming of this approach lies in the uncertainty of the exact distribution of scattering material along the 
line of sight (Lazio & Cordes 1998; Yusef-Zadeh et al. 1994), which in turn affects the amplitude and time scale of 
the predicted variations. To encompass the range of variations possible we consider a model in which the scattering 
material is entirely located in a single thin screen, either at a distance of 50 or 500 pc from Sgr A*, and one in which 
the material is distributed in an extended medium near Sgr A*, again with a scale length of either Az = 50 or 500 pc. 

Only variations caused by refractive interstellar scintillation are investigated here. Fluctuations caused by diffractive 
scintillation are possible in principle, but the extremely strong scattering observed toward Sgr A* renders the time scale 
of such scintillation of order seconds at centimeter wavelengths. Such variability is expected to be strongly quenched, 
given that recent estimates suggest the intrinsic source size far exceeds the angular scales probed by this phenomenon 
(Bower et al. 2004). 


3.1. The scattering model 
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The distribution of scattering material along the line of sight to Sgr A* is described in terms of the power spectrum 
of electron density fluctuations. This is modeled in the following standard form: 


/ 2 \-/3/2 

^iVe(q, z) = C%{z) + Rqlj exp 



( 2 ) 


where the amplitude of the power spectrum is written as a function of the spatial wavenumber, q, and the distance, 
z, from the source. The quantity Iq is the inner scale of the fluctuation spectrum, and is usually identified with the 
turbulent dissipation scale. 

The C 2 axial ratio observed in the scatter-broadened image of Sgr A* (Lo et al. 1998) indicates that the amplitude 
of the power spectrum varies with direction on the sky, presumably reflecting the orientation of the local magnetic 
field. The observed anisotropy is the result of either a change of the inner scale or the amplitude of the power spectrum 
as a function of orientation on the sky. We concentrate on the latter case because such anisotropy is expected for 
MHD turbulence (Goldreich & Sridhar 1995). In eq. l(21) the anisotropy is characterized by means of the parameter R, 
and is oriented so that the major axis of the scattering disk is along the x-axis. The anisotropy of a scattered image, 
C, is equal to the anisotropy parameter R when the length scales probed by angular broadening are much larger than 
the turbulent dissipation scale. On smaller scales the two measures are approximately equivalent (see Hd.2l belowb 

The spectral index of the turbulence is assumed to be independent of z, but its amplitude is allowed to vary through 
the quantity C%{z). For a thin screen located at a distance zq from Sgr A* one writes. 


rC^,|^-zo| < AL/2 
\ 0, otherwise. 


( 3 ) 


where the thickness of the medium, AL, is assumed to be far less than the source-observer distance, L. In the thin 
screen model the scattering measure SM = C^AL and the screen distance are the free parameters. 

When the scattering measure is large, as it is in the Galactic Center environment (see Hd.2l belowb the refractive 
modulations from a thin screen can be small. To this end, we also investigate a model in which the turbulent fluctuations 
are extended along the line of sight. Extended media are capable of producing larger refractive modulations relative 
to thin screens (Romani, Narayan & Blandford 1986; Coles et al. 1987). We consider the specific distribution in which 
the amplitude of the turbulence declines as a Gaussian from the Galactic Center: 

Clf{z) = (7^(0) exp ■ (4) 


The free parameters of this model are the scale length of the distribution, Az, and a normalization constant that 
sets the overall amplitude of the turbulence, (7^(0). The effective scattering measure of the medium in this model is 
SM = (7^(0)7r^/^Az/2 when Az is much smaller than the distance to the source. 

Unfortunately, the additional complexity inherent to the treatment of fluctuations in an extended medium forces us 
to abandon the explicit inclusion of anisotropy for this model. We take i? = 1 and normalize the model to a scattering 
strength intermediate to that implied by the scatter-broadening along the major and minor axis of Sgr A*. We feel 
that our failure to take anisotropy into account in a self-consistent manner in this model is a minor point compared 
to other uncertainties, such as the true distribution and strength of scattering material along the line of sight. In 
any case, the correct manner in which to incorporate anisotropy in an extended scattering medium model is highly 
uncertain. Goldreich & Sridhar (1995) point out that that the local value of i? in a thin plane of MHD turbulence, 
where “thin” is less than the outer scale of magnetic field fluctuations, is expected to be > 10^, and the value of 
anisotropic image broadening that one actually observes is much lower only because it represents an average over 
many different orientations of the magnetic held along the line of sight^. 


3.2. Determination of appropriate to the models 

The angular size of the scatter-broadened image of Sgr A* provides a direct means to determine the scattering 
strength appropriate for our scattering models. For a source of unit intensity whose intrinsic angular size is much 
smaller than the angular broadening size, the visibility^ observed on a baseline s is 


r(s) = exp 


1 

2 


D'^{sz/L,z) dz 


( 5 ) 


The quantity D'^ is the derivative of the phase structure function with respect to z and is given by 

/ OO 

^We(q:2:)(l-e*‘’'^)d2q. 

-OO 


( 6 ) 


^ The Goldreich ^ Sridhar (1995) theory of Kolmogorov MHD turbulence suggests that the local anisotropic ratio is 
where Lout is the outer scale of the turbulent magnetic field fluctuations and is the wavenumber on which the power spectrum is probed. 
^ We denote this observed visibility by P to distinguish it from the intrinsic source visibility, V used below. 
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The baseline Sg at which the visibility falls to 1/e of its maximum value at wavenumber k = 27r/A is related by 
9o = 2/kso to the angular radius at which the brightness distribution falls to 1/e of its maximum value, Oq. One 
therefore solves D'^{2z/k6oL, z)dz = 2 to determine the scattering strength. 

Measurements of angular broadening probe the structure function on a length scale ^ sz/L = 2z/k6QL^ which is of 
order kilometers for the case of Sgr A*. Consider, for instance, the scale probed by angular-broadening at 6 cm, where 
the angular diameter of Sgr A* is 49.6 mas along its major axis. Assuming the scattering occurs at z = 100 pc and 
taking the distance to the Galactic center as L = 8.5 kpc, one finds that angular broadening is sensitive to structure 
on scales of only 2 x 10^ m. 

This scale is much smaller than the expected inner dissipation scale. Spangler & Gwinn (1990) argue that the inner 
scale of the turbulent cascade is plausibly identified with the larger of the ion inertial length 


or the ion Larmor radius 



km. 


"^th „„„ 


T,; 


103 K 


1/2 



(7) 

( 8 ) 


where Ti is the ion temperature. The dissipation scale is larger than several kilometers for the range of plausible 
densities, temperatures and magnetic fields in the Galactic Center. 

The length scale on which angular broadening probes the scattering medium is important in determining the scatter¬ 
ing strength because the character of the phase structure function changes when its argument falls below the turbulent 
dissipation scale. The structure function scales as above this point and as below it. Defining Sdiff as the length 
scale over which the rms phase change is one radian, one sees that Sdiff oc 9^^. One has Sdiff oc for Sdiff > lo 

and Sdiff oc for Sdiff < h- Centimeter-wavelength observations show that the apparent size of Sgr A* scales as 
A^ oi±o. 03 y > 2 cm (e.g. Lo et al. 1998, Bower et al. 2004), which indicates either /3 « 4 or sz/L < Iq. However, in 
view of the foregoing arguments, only the latter explanation is viable. Scatter-broadening measurements at centimeter 
wavelengths yield no information on the spectral index of the electron density fluctuations. 

We normalize each scattering model by finding the appropriate Cf^AL or (7^(0) to reproduce the observed angular 
radius of the image of Sgr A*. In the thin screen model the phase structure function at r < Iq is approximated by 


D^{y)k - 
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2 ’ 2 ’ 2 
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(9) 


This equation is solved in conjunction with eq. © for a given R, z, L and 9q to find the scattering measure required 
to reproduce the angular size and ellipticity of the image of Sgr A*. In the regime r Iq applicable to the scatter¬ 
broadening of Sgr A*, the value of R and the image ellipticity do not correspond exactly, so one hnds the value of R 
numerically. 

In the extended medium model the integral D{sqz /L, z)dz can be inverted when Az < L and sqAz/L < Iq to 

find Cff{0) in closed form: 


C'w(O) 


V^r2A"‘Az3r(l — a/2) 


( 10 ) 


The models described below are normalized using the above expressions so that the scattering strength is always 
consistent with the observed degree of angular broadening toward Sgr A*. 

Although angular broadening measurements leave the spectral index of turbulent electron density fluctuations un¬ 
constrained, we concentrate on only two specific cases that illustrate the range of behavior possible, /3 = 11/3 and 
/3 = 3.9. The former case corresponds to Kolmogorov turbulence, while the latter case is of interest because steeper 
spectra produce refractive modulations whose amplitude exhibits a weaker frequency dependence. Such an index is 
suggested by the data, as the measured structure functions (Fig.© all saturate at similar amplitudes. Qualitatively 
different refractive scintillation behavior is possible for yet steeper, /3 > 4 power spectra (e.g. Blandford & Narayan 
1985, Goodman & Narayan 1985), but we do not consider such spectra here. Scintillations in this regime cause a large 
degree of image wander which is incompatible with observed limits placed on this effect by VLBI astrometry on Sgr 
A* (Backer & Sramek 1999; Reid et al. 1999; Reid & Brunthaler 2004). 
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Fig. 6. — Intensity structure functions from a thin-screen of scattering material with a turbulent power spectrum with index [3 = 11/3 
or 0 = 3.9, located a distance of 50 pc or 500 pc from Sgr A*. The strength of the scattering is set by requiring the scattering to reproduce 
the observed degree of angular broadening of Sgr A*. Because the turbulence is anisotropic the form of the intensity structure function 
depends on the orientation of the anisotropy axis with the scintillation velocity. The two limiting cases are shown for each wavelength: the 
curve which reaches the highest amplitude results when the velocity is oriented para llel to the major axis of the scattering disk, while the 
lower curve results when the velocity is oriented along the minor axis. Equation should be used to convert the values on the abscissae 
to units of time. 


3.3. Refractive Variations from a Thin Screen 

We quantify the amplitude and time scale of the intensity fluctuations expected from a thin scattering screen using 
the intensity autocovariance function, Ci{t). This quantity is directly related to the intensity structure function, 
Di{t) = 2[(7/(0) — Ci{t)], so it permits direct comparison between the models and the observed variations. 

Since the scattering medium is anisotropic the amplitude of the intensity variations on a given time scale depends 
on the relative orientation between the anisotropy and the scintillation velocity. As the direction of the scintillation 
velocity is unknown, we calculate intensity structure functions both when the velocity is parallel and perpendicular to 
the anisotropy axis. These two choices correspond to vjss = {v, 0) and vjss = (0, v) respectively. 

The intensity autocovariance due to scattering from a thin screen is (Codona et al. 1987, Coles et al. 1987), 


(7/(T) = 8^r2A2SM J d\(^^+Rql 


-/3/2 
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. z T viss • q 


where the phase structure function for an anisotropic screen is given by eq. for r <C ^0 and 

D^{r) = ^2 ^2 SM r(^l^-^/3/2) ^ ^2^2)/3/2-l^ 2 < /3 < 4, 


( 11 ) 


( 12 ) 


in the regime r > Iq. The function V(r) is the source visibility. We assume that the source is point-like and henceforth 
take V = 1. In practical terms, the visibility of the source is only important consideration in the thin-screen model 
if the intrinsic source angular diameter exceeds the angular broadening size. The assumption of a point-like source 
is valid at centimeter wavelengths, as Lo et al. (1998) and Bower et al. (2004) report that intrinsic source size only 
becomes comparable to the scatter-broadening size at wavelengths shorter than 7 mm. 

The scintillation velocity vjss depends on the motions of the source and scattering material relative to the observer, 
Vgrc.o and Vscreen.o respectively by (Gupta, Rickett & Lyne 1994): 



viss = 


(13) 
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P 

distance (pc) 

3.67 

3.9 

50 

SM= 2.70 X 10^5 

SM= 6.35 X 10^^ m-5-« 

500 

SM= 2.70 X 10^3 jjj-5.67 

SM= 6.35 X 10=1 m-®-9 



P 

screen thickness (pc) 

3.67 

3.9 

50 

C^(0) = 1.10 X 10'^ 

C^(0) = 2.52 X 10^ 

500 

C^(0) = 1.10 X 10'‘m-®'6'^ 

C^(0) = 2.52 X lO^m-®-® 


TABLE 2 

The scattering parameters used for each model investigated. All models assume a distance of L = 8.5 kpc to Sgr A* and 

AN INNER SGALE Iq = 100 KM. In THE THIN SCREEN MODEL THE SCATTERING PARAMETERS ARE DERIVED USING THE OBSERVED 
scatter-broadened DIAMETER OF SGR A* AT A6.03GM OF 49.6 MAS ALONG THE MAJOR AXIS. In THE EXTENDED MEDIUM MODEL, WHICH 
DOES NOT TAKE SCATTERING ANISOTROPY INTO ACCOUNT, THE SCATTERING PARAMETERS ARE MATGHED TO THE AVERAGE DIAMETER, 

37.5 MAS. 


Thus we see that the speed with which the scintillation pattern traverses the Earth is a factor ~ L/z larger than the 
speed of the scattering material itself. This correction factor is in the range 17— 170 for the plausible range of screen 
distances from Sgr A*. 

Numerical integration is used to derive the structure function given the scattering measure, anisotropy ratio - both 
derived from angular broadening measurements - and a screen distance and spectral index, (3. The results for both 
(3 = 11/3 and (3 = 3.9 and z = 50pc and z = 500pc are shown in FigureEl Table |21 lists the various models used and 
their associated parameters. 


3.4. Refractive Variations from an Extended Medium 

In an extended scattering medium the autocovariance of intensity fluctuations is given by (Codona et al. 1987, Coles 
et al. 1987) 


nL nCa 

Ci{t) = 87rrlX'^ / dz d^q 
Jo J — oo 
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/q(L-^ 


(q, z) sin^ 


q^z(L - z) 


2kL 
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,t] dt + i 


.^q-vissr 


(14) 


where 


h{t, z) 


f t{z/L-l),t < z, 
\z{t/L — l),t > z. 


(15) 


The assumption that the source is point-like is somewhat more problematic when the scattering medium extends 
close to the source, since clearly here the apparent angular source size must be appreciable relative to the local size of 
the scattering disk. Moreover, it is possible for the source size to affect the amplitude of intensity fluctuations without 
exerting a strong influence on the shape of the scatter-broadened image. This is because scattering material closer to 
the observer influences the angular broadening more strongly; eq. shows that one requires a smaller baseline s for 
a larger value of z to satisfy szjL = Sdiff- 

The effect of source size is investigated here by calculating the structure functions for a variety of intrinsic source 
sizes, from 1 /ras to 10 mas. The source is modeled as a circular disk whose visibility falls to zero on baselines longer 
than 1 /kOsrc so that the term V (q(T — z)/k) appearing in the integrand above is approximated as H[1/{L — z)0src — ?], 
where H is the unit step function. We also make the simplifying assumption that the scattering is always strong so 
that, even for very small z, the exponential in eq. d cuts off the integrand before the sine-squared term begins to 
oscillate (i.e. rdiff(z) < sfzjki). We therefore expand the sine function for small argument and obtain the following 
expression for the intensity autocovariance from an extended scattering medium 


= dzz^{L-zfe-^"l^^" dqq^-(^ 
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_{L - z)9s 


X exp 
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V k ■ 


,t] dt 


(16) 


The integrand over q cuts-off sharply once the source size becomes important, at (L — z) ^0sTc^ O'" when either of the 
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Fig. 7. — Theoretical structure functions as a function of wavelength for the extended medium model for a Kolmogorov spectrum of 
turbulence (/3 = 11/3) and for medium depths of 50 and 500 pc. The effect of finite size is illustrated by the various curves in each figure: 
red, brown, dark green, light green and blue dots represent the structure functions for a source of intrinsic angular size 10,1,0.1,0.01 and 
0.001 mas respectively. A finite turbulent dissipation scale of 10^ m is assumed in the calculation. 


two arguments of the exponential exceed unity, at 21^ ^, or which is defined by the implicit equation 
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where we write C = L/Az. Noting that the scintillation velocity is well-approximated by |vsrcreen,o — '^stc,o\L/z for the 
values oi Ljz ^ 1 under consideration here, the argument of the Bessel function in eq. lltill becomes qVgsAt, where 
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Fig. 8.— The same as FigiirelTihiit, a steeper, /3 = 3.9 electron density fluctuation power spectrum. 


^eff = |vscreen,o — Vsrc,o|- The intensity autocovariance thus reduces to 
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(18) 


Below we also use eq. CHI) to predict the amplitude of fluctuations at millimeter wavelengths. Since the scattering 
is weaker at shorter wavelengths, an additional wavenumber cut-off, at the inverse of the Fresnel scale, (^max = 
y/AirL/Xz{L — z), is introduced so that the amplitude of the variability is not overestimated. 

Equation Cl is integrated numerically to derive the intensity structure function expected due to refractive scintil¬ 
lation in an extended medium. These functions are shown in Figs. E|and|Sl 
The time scale at which the structure functions saturate can be understood in terms of the time required for the 
turbulent medium to traverse the scattering disk, of order L9q/2viss- A more rigorous estimate of the time scale is 
obtained by comparing the spatial wavenumber of the cut-off in eq. d with the form of eq. ©• The power spectrum 
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1.3cm 2cm 




Fig. 9.— A comparison between the observed structure functions and those predicted by the extended-medium A = 500 pc, (3 = 3.9 
model for a plausible range of source sizes, and the (3 = 3.9, 2 = 500 thin-screen model. A scintillation speed ujsg = 1000 km s~^ is assumed 
in this plot. 


of intensity fluctuations cuts-off when the term D'^(qh(t, z)lk^ t)dt reaches unity. The smallest wavenumber cut-off 
occurs at the outer boundary of the scattering medium, when z = AL and h{t,z) « t. Eq. 0, which describes 
angular broadening due to the scattering medium, involves an integral of similar form. By equating arguments in the 
two integrals one sees that the power spectrum cuts-off at g = 2/9oL, where do is the scatter-broadened size. This 
corresponds to a time scale r = Ldo/2uscreen- 

The intensity variations potentially yield information on the source size on scales well below the size of the scatter- 
broadened image. For instance, observations at 1 cm could distinguish between a source size of 1 /ias and 10 ^as simply 
on the basis of its variability properties. This might seem surprising because such small angular diameters represent 
only 0.1 and 1% respectively of the angular diameter of the scattering disk at this wavelength. This sensitivity to detail 
well below the scatter-broadened size is possible because phase fluctuations along the path of propagation are weighted 
differently between observations of scatter-broadening, which measures a second-order moment of the wavefield, and 
intensity variability, which represents a fourth-order moment. 

The ability to distinguish between such small source sizes arises in any scattering medium that extends very close to 
the source. Refractive modulations are strongest when the scattering strength is near unity. The scattering strength 
can be expressed as the ratio of the Fresnel scale rp = •\/z{l — z/L)lk to the transverse scale in the scattering medium 
over which the phase changes by one radian (i.e. the length scale, rdis, at which D^{rdis, z) = 1). The scattering 
strength increases with distance from the source as both the Fresnel scale and the amount of scattering material 
encountered along the ray path increase. Thus material close to the source can contribute greatly to the intensity 
variations provided that the source is sufficiently small that it does not substantially quench this contribution. 

4. THE DISTINCTION BETWEEN SCINTILLATION-INDUCED AND INTRINSIC VARIATIONS 

We now consider to what extent scintillation can account for the observed variability. To qualify as a viable 
explanation of the variability on any given time scale it must account for a large fraction of the modulation amplitude 
at that time scale. 

By this criterion no thin-screen model constitutes a viable explanation of the intensity variations in Sgr A* on any 
time scale. Larger variability amplitudes are predicted in this model when the scattering material is placed further 
from the source, but even when the screen is placed 500 pc from Sgr A* the variability amplitude is more than a factor 
of two below that observed. 

Another shortcoming of the thin screen model is its failure to account for the slopes of the observed structure 
functions at small time lags. None of the observed structure functions rise more steeply than r^, whereas the models 
predict that they should rise as Moreover, even if a thin-screen model were to reproduce the amplitudes and 

time scales at which the structure functions saturate it would still not explain any variations on shorter time scales. 
There does not appear to be any range of time scales for which any of the structure functions rise as steeply as the 
thin-screen model predicts. We thus conclude that thin-screen models constitute a poor explanation for the intensity 
variations observed on any time scale. 

Extended medium models fare significantly better at explaining the amplitude of the observed intensity variations. 
As can be seen from Fig. 13 models with shallow electron density power spectra {(3 < 11/3) still fail to account for 
most of the intensity variations, particularly at long wavelengths. However, the steep /3 = 3.9 power spectrum model 
reproduces both the saturation amplitudes ~ 0.1 of the observed structure functions and their weak dependence on 
wavelength. 

These structure functions also saturate at roughly the time scale observed in the data. The saturation time scales 
in the models can be understood as the time required by the scattering medium to traverse the spatial extent of the 
scattering disk. This rule of thumb estimates the time scale correct to within a factor of two relative to the time scales 
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indicated by the structure functions in Figs. [TjandlHl The predicted time scale is 

Certain extended-medium models also reproduce the generally shallow slope of the observed structure functions at 
small time lags. This behavior depends on the source size relative to the angular scales probed by the scintillation. 
Shallow slopes are present in some model structure functions over a large range of time scales, spanning up to two 
orders of magnitude in time scale before the structure function saturates. 

To illustrate the distinction between thin and extended-medium models, Fig.^ plots the observed 1.3 and 2 cm 
structure functions against several models. It is interesting to note that the (3 = 3.9 extended-medium model with 
an intrinsic source size of ~ 300 /ras would appear to closely match the observed structure functions. This size is 
comparable to that recently deduced by Bower et al. (2004) using VLBI. 

Our conclusion is that extended medium models with steep P > 3.9 power spectra are capable of reproducing the 
gross features of the variability at centimeter wavelengths. They do not explain all the detailed features of the structure 
functions. In particular, they fail to account for any of the various peaks apparent in the structure functions in the 
range 50 — 100 days which, as discussed in a which are probably due to flaring. Although it is hard to see how 
these features could be reproduced by a simple scintillation model, it is nonetheless pertinent to consider how the 
assumptions used in the models bear on the predicted scintillation properties. 

In all extended medium models it is assumed that all layers of the scattering medium move with identical peculiar 
velocities. When this is not the case different layers of the medium may cause intensity fluctuations on different time 
scales. The extent to which this could occur depends on the relative contributions that layers at different distances 
make to the intensity fluctuations. Significant variations will only be observed on different time scales if two layers 
which both contribute substantially to the scattering each possess different transverse velocities. It is possible to see 
which layers contribute most to the scattering by comparing the amplitude of the structure functions for medium scale 
lengths of 50 and 500 pc. The intensity fluctuations are dominated by scattering layers within 50 pc from Sgr A* if 
the amplitude on a given time scale is identical in the two models. Conversely, when they differ, most of the intensity 
fluctuations originate in the scattering medium beyond 50 pc from Sgr A*. 

The time scale may also vary from the wavelength scaling predicted if different layers of the medium move at 
different velocities. This is possible because the layer that contributes most to the intensity fluctuations changes with 
wavelength (see Fig. EJ. 

Another assumption inherent to the scattering model lies in the simplicity of the source structure. Our models assume 
the simplest possible structure: a source comprised of one single circularly-symmetric component. More complicated 
structure could cause qualitatively different variability. To be specific, eqs. jnj and m establish how the power 
spectrum of the intrinsic source brightness distribution, |R(r)p alters the power spectrum of scintillation-induced 
intensity fluctuations. 

Consider, for instance, how a source comprised of two compact components (e.g. a jet and a counterjet) would alter 
the scintillation characteristics. Whereas the visibility of a single point source is constant and independent of baseline 
length, the visibility amplitude of a double source with angular separation A9 pointing along the direction 9 ~ A6/A6 
oscillates on a scale of length rose = 9/kAO. This oscillation enhances the scintillation fluctuations at certain time 
scales relative to others. The first peak of the oscillation emphasizes the power spectrum of intensity fluctuations 
at a fundamental spatial wavenumber q = 9[A0{L — z)]~^, for a scattering layer a distance z from Sgr A*. This 
corresponds to fluctuations on time scales tq « A6{L — z)/wscreen- The amplitudes of visibility oscillations at higher 
harmonics depend on the size of the components constituting the double source. The visibility amplitude of a double 
source comprised of sufficiently compact components contains peaks comparable to the amplitude of the fundamental 
peak at integer multiples of rose- This in turn would enhance the scintillation fluctuations at shorter time scales, with 
the nth harmonic enhancing fluctuations on a time scale t ro/n relative those on surrounding time scales. 

It is possible in principle for source structure to explain additional features of the fluctuations observed toward Sgr 
A*. One could appeal to a double source structure to explain peaks in the observed structure functi ons. However, 
this model cannot not explain the peaks at time lags between 30 and 150 days that are discussed in io A double 
source would give rise to a number of regularly-spaced peaks in the structure function, not a single, isolated peak. 
Moreover, the source separation required to explain the location of the peaks is sufficiently large that it would have 
been observed. A peak on a time scale of r = 50 days would require a double source of separation 3 mas for a scattering 
medium moving at 1000km s“^. 

We conclude that the interpretation of the isolated peaks present in the structure functions discussed in 42.21 a,s flares 
intrinsic to Sgr A* is robust to the assumptions made in the scintillation models considered here. 

4.1. The predicted role of scintillation at millimeter wavelengths 

Given the success of the /3 = 3.9 extended medium model in explaining the broad characteristics of the centimeter 
variability, we have applied it to the predict the variations at millimeter wavelengths. Sgr A* flux monitoring is already 
being carried out at the Sub-Millimeter Array (Zhao et al. 2003), and similar monitoring will soon be possible using 
CARMA and, eventually, ALMA. 

The predictions at 1 and 3 mm from the AL = 500 pc, /3 = 3.9 extended-medium model are shown in Fig. 

This model implies that a 1 /ias (10 /ras) source should exhibit 25% (22%) root-mean-square fluctuations on a time 
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Fig. 10.— The structure function of the predicted variations for the best-fitting, AL = 500 pc, /3 = 3.9, scintillation model at centimeter 
wavelengths applied at 1 and 3 mm. 


scale of 1.4 (Psjcreen/lOOO kms~^) hours at 3mm. A 100 ^as would exhibit only 2.6% variations, and on a time scale 
approximately three times longer. Given recent measurements of the intrinsic size of Sgr A* at 7mm of ~ 0.22mas 
(Bower et al. 2004) and assuming a size-dependence it is reasonable to expect scintillation variations of order 3% 
at a wavelength of 3 mm. 

The scintillation characteristics are even more sensitive to source size at 1 mm. Only a source size of ~ 1 ^as is 
sufficiently small to exhibit 26% fluctuations. Sources of 10 and 100 p,as would exhibit r.m.s. fluctuations of 3.3 and 
0.29% respectively. Again, the predicted variations occur on intra-day time scales of 4 and 12 hours respectively 
assuming a scintillation velocity of 1000 km s“^. 

It is possible to already compare these predictions with a number of observations. Both Wright & Backer (1993) 
and Tsuboi et al. (1999) report variations of order IJy amplitude at a wavelength of 3 mm. Such large flux density 
excursions are difficult to explain in terms of the present scintillation model, suggesting that intrinsic source activity 
is responsible for most of the variability. 

Mauerhan et al. (2005) have recently claimed the Sgr A* undergoes only ~ 20% intra-day variations at 3 mm. 
These variations, if real, would require an intrinsic source size of 10 — 30 ^as at this wavelength to be consistent with 
scintillation. 

5. DISCUSSION & CONCLUSIONS 

Our analysis of multi-frequency monitoring data presented by Zhao et al. (1992, 2001), Herrnstein et al. (2004) and 
Falcke (1999) indicates that Sgr A* exhibits no quasi-periodic oscillatory behavior on any time scale between one week 
and 200 days. The variability amplitudes are remarkably constant with frequency, varying between 30 and 39%, but 
the time scale on which they saturate increases with wavelength. 

Several structure functions show evidence for variability on multiple time scales. If the errors associated with the 
flux density measurements are correct, the data of Herrnstein et al. (2004) indicate that the source exhibits unresolved 
6-10% inter-day (< 4-day) variations between 7 mm and 2 cm. No structure functions exhibit, within the errors, any 
evidence for appreciable variability with time scales longer than 1000 days. The long term stability of the radio flux 
implies there is very little long-term variation in the accretion rate. 

We sought to explain the general features of the variability by reproducing the shape and amplitudes of the observed 
structure functions using several scintillation models. Both thin-screen and extended-medium models were considered. 
No thin screen model accounts for the properties of the variations. The structure functions of the observed lightcurves 
rise less steeply with time lag and saturate at higher amplitudes than predicted. They underpredict the amplitude of 
the variability by at least a factor of two at the saturation time scale and often by a more than an order of magnitude 
on shorter time scales. If the medium responsible for the scattering of Sgr A* lies on a thin screen all of the observed 
flux variability must be intrinsic to the source itself. 

Certain extended-medium models, on the other hand, do explain the amplitude of the fluctuations over a large range 
of time scales. Models in which the electron density fluctuations follow a Kolmogorov power spectrum, corresponding 
to P = 11/3, and a slightly steeper, /3 = 3.9, power spectrum were investigated. Of the two models, only those with 
a /3 = 3.9 spectrum account for the amplitude of the fluctuations at all wavelengths. If scintillation is to precisely 
predict the amplitude of the flux density variability of Sgr A*, this suggests that the power spectrum of the Galactic 
Center turbulence is slightly less steep, with an index lying in the range 3.8 < /3 < 3.9. 

The most successful extended-medium model examined was used to predict the maximum contribution that scin¬ 
tillation could make to future observations of Sgr A* at millimeter wavelengths. The expected variability amplitude 
depends strongly on the intrinsic source size. A 1 — 10 /ias object at 3 mm would undergo fractional root-mean-square 
fluctuations of ~ 25%, but a 100 /xas source would exhibit only 3% variations. Scintillation is even more sensitive to 
source size at 1 mm, with a 1 p,as source expected to display 26% variations, but a 10 /ras source would display only 
3% variability. 

With what physical structures can we associate the extended scattering medium? Given only the approximate 
match between our model structure functions and the actual structure functions, we do not expect that the extended 
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scattering medium must in fact be parameterized as discussed in the text. In fact, the extended scattering medium 
might consist of only a few thin scattering media distributed over a range of distances from Sgr A*. These thin media 
would have characteristics similar to the thin medium discussed by Lazio & Cordes (1998), with densities > 10^ cm“^. 
The details of the structure functions are not sufficient to constrain this result further. For comparison purposes, we 
do compute the mean density associated with the extended scattering medium as presented. In this case, the mean 
density is ~ 10cm“^, assuming the outer scale of the turbulence to be Ipc. This density is substantially greater than 
the density of the diffuse hot ionized gas (Tg ~ 10^ K, Ue ~ 0.05 cm“^) detected in the Galactic Center. As Lazio & 
Cordes (1998) discuss, the scattering medium may be the interface between this hot medium and molecular clouds in 
the central 100 pc. Our results are consistent with the picture reached by Lazio & Cordes (1998) for the scattering 
medium with the modification that the scattering may take place at a range of distances from Sgr A*. 
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